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In the framework of the time dependent local density approximation, applied to valence 
electrons, coupled non-adiabatically to molecular dynamics of ions, the irradiations of 
ethylene by laser and fast charged projectiles are studied. We find that the Coulomb 
fragmentation sensitively depends on the laser frequency and on the charge of the pro- 
jectile. 

1. Introduction 

It is well know that the interaction of ionizing radiation with biological tissue can 
induce severe damage to DNAP^ With the development of ion sources such as the 
electron cyclotron resonance (ECR) ion source and modern accelerators such as cool- 
ing storage ring (CSR), a lot of investigations have been devoted to understanding 
the molecular mechanisms underlying biological radiation damage. More recently, 
a number of experimental investigations have been addressed concerning the action 
of ions on molecules with high and low energies From the mass spectra, the 
relative cross sections for the different ionization and fragmentation channels can 
be evaluated. One can observe not only the process of electron capture but also the 
highly excited processes, such as electronic emission and Coulomb explosion due to 
the induced ionization. 

Impact of radiation also includes excitation by photons and we will discuss 
their effects as well. The advantage of photons is that modern lasers provide flex- 
ible and powerful pulses for studies of photo-induced dynamics, such as optical 
response,^ multi-photon ionization,^ and Coulomb explosion of clusters,^ high- 
order harmonic generationp^Dbond softenmgj 19 * 20 land charge resonance enhanced 
ionization! 21 1 22 1 
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The perfect theoretical description of these mechanisms would require to solve 
the time-dependent Schrodinger equation for all electrons and all nuclear degrees of 
freedom. However, this ultimate attack exists only for very small systems, such as 
atoms p2EHH2 and H J !2SHH] F or larger systems, more approximate ways have to be 
found. Recently, Horsfield et al. used correlated electron-ion dynamics to study the 
excitation of atomic motion by energetic electrons.^ Saalmann et al. have devel- 
oped a non-adiabatic quantum molecular dynamics (NA-QMD) to study different 
non-adiabatic processes in different systemJSKsH and Calvo et al. used a combined 
method to study the fragmentation of rare-gas clusters.^ It is also interesting to 
develop a fully-fledged coupled ionic and electronic dynamics. At the side of elec- 
tronic dynamics, the well tested Time-Dependent Density Functional Theory at the 
level of the Time-Dependent Local Density Approximation (TDLDAJPSl i s used. The 
ions are treated by classical molecular dynamics (MD). This altogether provides a 
coupled TDLDA-MD, for a review seeP^Up to now, this approach has been applied 
to study the dynamical scenarios of simple metal and hydrogen clusters f22HH] g e _ 
sides these systems, organic molecules are particularly interesting cases and recently 
motivated a lot of theoretical works to investigate cither the optical response in the 
linear domairP^ or the dynamics of it without and with considering the external 
laser fieldPEH 

In this paper, the irradiations of ethylene by laser and fast charged projectiles 
are studied by using a coupled TDLDA-MD 2 ^ and various non-adiabatic effects are 
investigated. The paper is organized as follows. Section 2 provides a short presen- 
tation of the theoretical and numerical approach. Section 3 first gives the optical 
response properties of ethylene and the effect of ionic motion on the excitation dy- 
namics is discussed. Then different excitation scenarios of ethylene as a function 
of the laser frequency and charged projectiles are presented. Section 4 gives some 
conclusions. 



2. Theory 

In this section, we briefly represent the real-time method in TDLDA-MD. The 
molecule is described as a system composed of valence electrons and ions. Our 
test case, the ethylene molecule C2H4, has 12 valence electrons and 6 ions. The 
interaction between ions and electrons is described by means of norm-conserving 
pseudopotentialspSl 

Valence electrons are treated by TDLDA, augmented with an average-density 
self- interaction correction (ADSIC)P^ They are represented by single-particle or- 
bitals c/>j(r,t) satisfying the time-dependent Kohn-Sham (TDKS) equations,^ 

d ~ / V 2 \ 

i— 4>j(r,t)= H Ks <p j (v,t)= +V eS (r,t)\<j> j (r,t), , j — 1, N. (1) 

The term V e g stands for the Kohn-Sham effective potential and is composed of four 
parts, 

V eS (r,t)) = V ion {r,t) +KxtM) + V n [n](r,t) +V KC [n](r,t), (2) 
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where V\ on — J2i ^s( r — R-j) is the ionic background potential, V cx t the external 
potential, Vh stands for the time-dependent Hartree part and V xc is the exchange- 
correlation (xc) potential. The electron density is given by n(r,t) — \<fij(r,t)\ 2 . 
The xc potential 14c [n] (r, t) is a functional of the time-dependent density and has 
to be approximated in practice. The simplest choice consists in the TDLDA, defined 
as 

X homf 

Vg D ^[„](,,t)= 



(3) 

=n{r,t) 

where £x° m [n] is the xc energy density of the homogeneous electron gas. For e^°" 



xc 



we use the parameterization of Perdew and WangPS The form of the pseudopo- 
tential for a covalent molecule is taken fromj^&l including a non-local part. The 
TDLDA approximation is augmented by an Average Density Self-Interaction Cor- 
rection (ADSIC) 47 to put the single-particle energies at their correct values. 

The ground state wave functions are determined by the damped gradient 
methodP The TDLDA equations are solved numerically by time-splitting tech- 
nique.^ The nonlocal part contained in V c fr is dealt in an additional propagator 
and treated with a third-order Taylor expansion of the exponential.^ Absorbing 
boundary conditions are employed to avoid reflecting electrons.^ 

Ions are treated classically and propagated by classical Molecular Dynamics 
(MD) equations 

m I ^±=F I {R I ,t) (4) 

The external perturbation by the laser field, neglecting the magnetic field com- 
ponent, acting both on electrons and ions, is given as 

Kxt.las = E zfi as (t) cos(tJi as £) (5) 

where Eq oc vJ, I denoting the intensity, z is the dipole operator, f\ as (t) is the 
pulse profile and wi as is the laser frequency. In this paper, /i as (i) is chosen as a cos 2 
in time. 

The external perturbation by a charged projectile is modeled as 

Qe 2 

>4xt,ion — , = (6) 



\J[x~ (v P rojt + x )] 2 +y 2 + (z - b) 2 



where Q, fproj, £o denote the charge, the velocity and the initial position of the 
projectile respectively; b is the impact parameter. The excitation amplitude may be 
tuned through variation of Q, w pro j, and b. Here we focus on the influence of Q, b 
and collision orientation. The velocity is fixed to u pr oj = 20a /fs. 

The calculations are carried out using a cubic box of size 72 x 72 x 64 with 
a grid spacing of 0.41 ao and a constant time step At = 6 x 1CP 4 fs. The dipole 
moment can be obtained by Di{t) = J d 3 r rj n(r, t), with i — x,y,z. The number 
of escaped electrons is defined as N csc = N(t = 0) — J v d 3 r n(r, t), where V is a 
sufficiently large volume surrounding molecule. A detailed link with experiments is 
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the probabilities P k (t) of finding the excited molecule in one of the possible charge 
states k to which they can ionize. The formula are deduced from the occupied 
numbers of the time-dependent single-particle orbitals and are taken from.^ 

3. Results and discussion 
3.1. Basic properties of C2H4 

Fig. [T] shows the spatial structure of ethylene. The molecule is planar and we place 



Fig. 1. Ionic structure of 




ethylene. 

it in the x-y plane with two carbon atoms on the x axis and four hydrogen atoms 
in x-y plane. The center of mass is at the origin. The single-electron energies range 
from —23.3 eV for the deepest bound valence level to the HOMO at —11.5 eV. The 
latter is in good agreement with the experimental ionization potential (IP) of 11 
eV^ (relative error of 4.5 %). Similar errors are obtained when comparing the C-C 
and the C-H bond lengths (2.43 and 1.97 do respectively, while the experimental 
values are 2.53 and 2.05 acpO). 

Fig. [5] shows the optical response, i.e. the photo- absorption strength distribu- 
tion, of C2H4. It shows strong, isolated peaks at lower energies and changes to a 
continuum above the IP (11.5 eV). This result is consistent with former TDLDA 
calculations.^ The lowest peak is at 7.74 eV, in good agreement with previous CI 
calculations (7.76 eV^). We obtain the highest peak in the optical response lying 
at 8.16 eV. The frequency selective laser pulses will be sensitive to this much fluc- 
tuating spectrum while a collision with fast ions will deliver a broad spectrum of 
frequencies and will excite all modes at once. 

3.2. Laser excitation 

3.2.1. The key role of non-adiabatic couplings 

To investigate the effect of the motion of the ionic centers of the molecule on the 
excitation dynamics, we compare, for the case of a laser excitation, calculations 
either with full ionic motion or with fixed ionic configuration. The laser intensity 
is 10 13 W/cm 2 and its full width at half maximum (FWHM) is 100 fs. We have 
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Fig. 2. Total 
photo-absorption strength 
of C2H4. The vertical ar- 
row corresponds to the 
value of the ionization po- 
tential of the molecule. 



chosen the laser frequency equal to 8.16 eV, that is the value of the highest peak 
in the optical response (see Fig. [5]). We thus expect a resonant excitation of the 
ethylene. Fig. [3] shows the time evolution of the x ionic coordinates, the number of 
escaped electrons (N osc ), and the dipole moment in the x direction, D x , for moving 
and fixed ions. Up to about 35 fs, no significant difference is visible in D x and N esc . 
From then on, comparing both cases, the dipole moments diverge rapidly in time 
and more complex pattern are observed for the case of moving ions. Consequently, 
the total ionization also differs and more electronic emission is obtained when ions 
can move (0.8 additional electron is extracted). The ions precisely start to exhibit a 
sizeable motion at about 40 fs (see top panel). The C-C bond oscillates (see the open 
squares) while the H (close squares) escape the system, due to the large charge. The 
non-adiabatic coupling between ions and electrons allows a feed-back of the ionic 
motion on the electronic response, as can be tracked on the dipole oscillations. The 
enhancement of ionization for the moving ions case is related to the expansion of the 
molecule. Indeed, after increasing ionization, the optical spectrum is blue-shifted, 
so aji as becomes smaller than the resonant frequency. However, the Coulomb force 
from increased net charge impels the expansion of ions which red-shifts of the optical 
spectrum and allows u;i as to be in the region of resonance again. The difference of 
electronic emission above shows the importance of treating electronic and ionic 
dynamics simultaneously. 

3.2.2. Scanning laser frequency 

To explore in more detail the influence of the laser frequency, we consider two cases, 
namely wiascr = 9.5 eV, lying is in the resonant region, and wi aS cr = 6.8 eV which 
is off-resonant. Both cases are depicted in Fig. |4]in the left and the right columns 
respectively. The laser parameters are a polarization along x axis, an intensity of 
I = 10 14 W/cm 2 , and a rather short FWHM of 20 fs. The bond lengths and the 
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Fig. 3. Excitation of ethylene in- 
duced by a laser pulse (I = 
10 13 W/cm 2 , uj laB = 8.16 eV and 
FWHM of 100 fs), for fixed or 
moving ions. Top panel: time evo- 
lution of the x ionic coordinates, 
when ionic motion is allowed. Mid- 
dle and bottom: iV eS c and dipole 
moment D x as a function of time, 
for moving (dark lines) or fixed 
(green curves) ions. 
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r.m.s. radii in the resonant case, top-left panels of Fig. HI show much different 
dynamics for the various bonds. The C-H bonds grow quickly and monotonously, 
which means that the bonds are broken and hydrogen ions are flying apart in the 
x — y plane. The C-C bond (dark solid line), on the contrary, oscillates with a rather 
large amplitude during the whole time span. The off-resonant case, top-right panels 
of Fig. |4j shows regular and faint oscillations for all bonds and r.m.s radii. The 
molecule remains unfragmented, although highly excited. 

The ionization N csc and the dipole evolution shown in the bottom panels of 
Fig. HI corroborate these findings. In the off- resonant case (right panels), the dipole 
nicely follows the laser field, producing an emission of 0.6-0.7 electron during the 
pulse duration. As soon as the laser is switched off, faint dipole oscillations remain 
and Nesc levels off. The resonant case visibly produces a complex pattern in the 
dipole evolution with large amplitudes; more ionization is observed, driving the 
molecule above the limits of binding. Note that -/V esc starts to increase at around 5 
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Fig. 4. 

Excitation of ethy- 
lene 

by a laser pulse of 
I = 10 14 W/cm 2 , 
FWHM of 20 fs, 
polarization along 
x axis, for two dif- 
ferent frequencies, 
Wiao = 9-5 eV (left 
col- 
umn), and aJi as = 
6.8 eV (right col- 
umn) . From top 
to bottom panels: 
time evolution of 
bond lengths, ex- 
tensions of 
the ionic distribu- 
tion in three direc- 
tions, dipole signal 
along the laser po- 
larization and the 
number of escaped 
electrons. 



fs while the ionic motion (see the bond lengths) follows later. Ionization is directly 
caused by the laser pulse and electrons react quickly because they have a low mass. 
The ions react indirectly, namely to the Coulomb pressure induced by ionization, 
and more slowly due to their larger mass. 

The above discussion demonstrates the versatility of a laser excitation when 
playing with its frequency (at low intensities). Indeed the value of this frequency 
plays an important role in the excitation dynamics of ethylene. In the resonant 
region, enhanced ionization often leads to a Coulomb fragmentation of the molecule. 
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3.3. An example of ionic collision 

We finally end with a typical example of collision between ethylene and a projectile 
with charge Q = 2. The projectile moves along x direction at the velocity of 20 ao/fs, 
with an offset of b — 7 ao in z direction. The initial position of the projectile 
is (— 100 ao, 0, 7a ), thus safely far away from the ethylene. Fig. [5] displays the 
time evolution of the dipole moments in the three spatial directions, N csc and the 
ionization probabilities p( n +). Let us start with the dipole response after collision 



Fig. 5. Time evolution of the dipole 
moments (top), of the total number of 
escaped electrons (middle) and ioniza- 
tion probabilities (bottom), after col- 
lision of C2H4 with a projectile of 
charge Q = 2, impact parameter b = 
7ao and velocity along the x direction 




of v„ 
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with the projectile (top panel). Since the projectile moves in the x — z plane along 
the x axis, no significant excitation along the y direction is expected. And indeed, 
the dipole moment in that direction, D y , is vanishingly small. On the contrary, D x 
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and D z exhibit at time of closest approach (at about 5 fs) an almost instantaneous 
large shift and quickly relaxe towards much more gentle oscillations that last for a 
long time. Note that the shift in x direction is higher than that in z direction. This is 
consistent with the fact that the extension of the electronic cloud is larger in x than 
in z directions (see the ethylene configuration in Fig.[T]). The total ionization, shown 
in the middle panel of Fig. [5j evolves according to the dipole excitation : it jumps 
almost immediately at impact time in correlation to the huge dipole moments, 
while the faint later dipole oscillations do not induce further electron emission. 
The bottom panel of Fig. [5] shows the detailed ionization probabilities. Their time 
evolution proceed similarly to that of N csc . Indeed p( 0+ ) falls down abruptly at 
the time of closest approach, while the other charge state probabilities grow very 
quickly. Then all probabilities level off and no evolution in time is observed. The 
charge state 1 + has the largest probability in accordance with a net ionization of 
nearly 1 (see middle panel), but a sizeable value of P 2+ of 20 % is nevertheless 
attained. As for higher charge states, their probabilities are negligible as expected. 

At the side of the ions, this fast collision only induces a visible motion in the 
x — y plane. More precisely, faint oscillations parallel to the x axis of the C-C bond 
(aligned with the x axis) are observed, while the C-H bond exhibit larger (but still 
small) oscillations, because of the smaller mass of the H compared with that of the 
C. In the y direction, the C atoms simply do not move and the H gently oscillate 
around their equilibrium positions. 

Apart from the different dipole responses, this collision case is close to the off- 
resonant laser excitation presented in Fig. IH right column. Indeed, the total ion- 
ization there is quite similar (0.7 compared with 0.9 here). However, the ethylene 
molecule seems able to cope such a ionization and remains in one piece. 

Finally, we just mention here some other cases of projectile collision that we 
have explored. Increasing the charge of the projectile can cause a direct Coulomb 
fragmentation of the ethylene molecule. We have also varied the value and the 
direction of the impact parameter b. When the projectile moves perpendicular to 
the ethylene plane, the electronic response is of course weaker, since the interaction 
time is shorter than in the case of a projectile moving parallel to the molecular 
plane. Still the investifated test cases lead us to conclude that the dominant (first 
order) parameter is the actual (absolute) value of b. 

4. Conclusions 

In this paper, we have demonstrated the capability of TDLDA-MD to describe 
complex coupled ionic and electronic excitations. We discussed various excitation 
scenarios of ethylene subjected to different laser pulses and fast charged projectiles. 
These scenarios involve both electrons and ions but the relative role of each species 
does depend on the actual excitation conditions of the laser pulse and charged 
particles. The response in time of the electron cloud, such as dipole deformation, 
total number of escaped electrons and ionization probabilities have been presented. 
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Varying the laser frequency with fixed laser intensity and pulse length indicates that 
the appearance of de-excitation/explosion scenarios depends on the relationship of 
the laser frequency with the eigenfrequencies of the system. We have also explored 
the case of a projectile collision at moderate velocity, a situation that can be easily 
handled experimentally. A very fast excitation is reported. A systematic study of 
the dependence on charge, velocity and impact parameter of the projectile will be 
presented in a forthcoming publication. 
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